Effects of nonintegrability on stabilization of Feshbach molecules in atom waveguides 
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Deactivation of broad quasi-one-dimensional dibosonic molecules is analyzed. Within integrable 
Lieb-Liniger-McGuire (LLMG) model an exact expression does not demonstrate suppression of the 
deactivation at low collision energies. Solution of Faddeev equations demonstrates that when a 
Feshbach resonance lifts the symmetry of the LLMG model the deactivation becomes suppressed. 
This effect shows a way for formation of a stable gas of dibosonic Feshbach molecules. 
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Ultracold molecules have recently been formed by Fes- 
hbach resonance association (see reviews jPj and original 
works 0, [H ) . These Feshbach molecules are very broad 
since they are superpositions of the closed and open chan- 
nels. However, inelastic collisions of atoms and molecules 
lead to strong losses of molecules composed from Bose 
atoms in experiments Q ■ In the fermionic case 0] inelas- 
tic collisions are strongly suppressed due to Pauli block- 
ing and weak coupling of the atoms [3| ■ This effect allows 
investigation of BEC-BCS crossover and other phenom- 
ena of fundamental physical importance. 

Tight confinement of atomic motion in two directions 
by atom waveguides strongly modifies atomic collisions in 
quasi-one-dimensional (ID) regime, when the transverse 
excitation energy wj_ (in units where h = 1) substantially 
exceeds the collision energies [5|. This regime has been 
realized in 2D optical lattices @, 0|, elongated atomic 
traps , and atomic integrated optics devices [t| . Broad 
quasi-lD molecules, predicted in [l(| and observed in Q, 
have binding energies less then ui±. Such diatoms can be 
described by two-channel ID model 11, [Hj], where the 



closed channel incorporates both 3D closed channel and 
excited transverse waveguide modes. 

Rates of deactivation into tightly bound (non- 
Feshbach) states are approximately proportional to the 
probability to find three atoms (two of which belong to 
the molecule) in the same place, e. g. to three-body 
(3B) correlations 0]. 3B correlations were analyzed in 
[13j within exactly-soluble Lieb-Liniger model [lj|. In 
this model ID bosons interact by zero-range potentials 
U a S (zj — zi), where Zj are atomic coordinates and the 
interaction strength U a is energy-independent. In the 
case of repulsive interactions (U a > 0), which does not 
bound the atoms, the correlations are suppressed at low 
collision energies and strong interactions 131 ] . Bound 
states can be formed in the case of attractive interactions 
(U a < 0) described by the McGuire solution 15]. The 
tree-atom correlations for atom and diatom with relative 
momentum po can be represented, using 3B wavefunction 
ipo (zi, z 2 , z 3 ) 0, as 
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FIG. 1: Schematic description of 2B channel potentials. 



where m is the atomic mass. This expression has a 
non-zero limit at low collision energies (po — > oo) or 
strong interactions (\U a \ — > oo). This behavior differs 
from the free-atom case 



13| . since bound atoms keep 
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non-vanishing imaginary momenta ±^mU a even when 
Po — > 0. Thus, inelastic collisions of ID broad molecules 
are cardinally different from collisions of two structure- 
less particles, which are suppressed both within the Lieb- 
Liniger-McGuire (LLMG) model and beyond it 
This difference reflects the fact that broad Feshbach 
molecules can not be considered as zero-range objects 
and their deactivation is not a 2B process. 

The present Letter demonstrates that deactivation of 
broad ID molecules becomes suppressed when integra- 
bility of LLMG model is lifted by a Feshbach resonance. 
Similar effects can be expected for other mechanisms of 
integrability lifting, e. g. due to virtual transverse mode 
excitation [18|, [l!|. In contrast with these effects, other 
processes, such as reflection and dissociation in atom- 
diatom collisions and three-atom association 17] , become 
allowed when the nonintegrability is lifted. 

Except of demonstration of a new observable effect of 
non-integrability, which is interesting for atomic, molec- 
ular, and statistical physics, as well as for quantum field 
theory, present results show a way for formation of a sta- 
ble gas of dibosonic Feshbach molecules. 

Consider multichannel collisions of ID Bose atoms de- 
scribed by the annihilation operators ^ a (z). The model 
includes several 2B channels (see Fig. [IJ. The Feshbach 
closed-channel state, described by the molecular annihi- 
lation operator ^ c (z), lies at the energy D c close to the 
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open channel threshold, which serves as the energy origin. 
A set of deactivation product channels {ci}, described by 
the molecular annihilation operators w<j (z) , lie at the 
energies D d far below the open channel threshold. The 
system can be described by the Hamiltonian 



H = J dzk>t( 



z) 



ip (zi, z 2 , z 3 ), which can be reduced to the Faddeev- 
Lovelace equation for the symmetric transition amplitude 

ffl^ 1 f 

X(p,Po) = 2Z(p,p ) + — — I dqZ(p,q)T 1D (k(q))X (q,p ) 
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+^Lh (*) + E (*) + (z) + V cd (z) + (z)] 1.(2) 
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Here the interaction VFesh (z) = gW c (z) (z) ^ a (z) de- 
scribes the Feshbach coupling, and 

Vad (z) = d ad ¥ a (z) ¥ a (z) (z) , Vcd (z) = d cd ¥ c (z) (z) 

(3) 

are couplings of the open and closed channels, respec- 
tively, to the product channels. The non-resonant in- 
teraction strength £/ Q , resonance detuning D c , and the 
Feshbach coupling strength g can be related to parame- 
ters of atomic collisions and the waveguide llj, [l2| [see 
Eq. (JT5J below]. 
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Like in the two-channel case [17], 2B problem can be 
described by the ID T matrix 



Tid (k) 



Uj (k 2 /m) + -mk- 1 



which depends on the relative momentum k of two col- 
liding atoms, where the energy-dependent interaction 
strength U e s(E c ) incorporates effects of all channels. 
The poles of Tid (k) on the positive imaginary axis, 
k = in n , correspond to 2B bound states (diatoms) with 
energies — ft„/m. They are superpositions of the open, 
closed, and deactivation product channels. The diatoms 
have finite size (~ k^ 1 ), although this model approxi- 
mates the closed-channel and the deactivation-product 
molecules to be infinitesimal in size. 

Substitution of the state vector for the three-atom sys- 
tem 



Here k (q) = \JmE + iO — 3g 2 /4, E is the total energy in 
the center-of-mass system, p n = 2y (mE + k 2 ) /3 is the 
relative momentum of the atom and diatom in the state 
n, and n = corresponds to the initial diatom state. The 
probabilities of inelastic reflection and transmission with 
the diatom transition to the state n ^ can be expressed 
as 
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where k < n n and W n is the contribution of the open 
channel into the diatom state n (see 12|, |l7|). The deac- 
tivation rate coefficient is 



3pp 
2m 



[Prof {0^n) + P tr (0 -> n)] . (8) 



In the case of high deactivation energies, 
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T\d (k) has poles at k = in d « iy/ m\D d \, corresponding 
to the deactivation products. Other poles are approxi- 
mately determined by the same cubic equation as in the 
two- channel model [13, 0] , 



-U a K 



mD c K + —m 2 D c U a — m 2 \g\ 2 = 0. (10) 



Contributions of the product channels into the corre- 
sponding weaker-bound diatoms (k„ <C n dl n ^ d) can 
be neglected. Whenever U a > or D c > 2\g\ 2 /U a , when 
Eq. (fl~0|) has a single real positive root, the deactivation 
rate coefficient can be approximately expressed as 



|* 3 > = J d 3 zcp Q ( Z1 ,Z2, *a) *i (zi) *J (z 2 ) *t ( Z3 ) 

+ E / d zdz m cpd (z, z m ) (z) *t (z m ) 1 1 vac) ,(5) 

where |vac) is the vacuum state, into stationary 
Schrodinger equation with the Hamiltonian |(5J) leads 
to coupled equations for the wavefunctions of the 
open (three-atom) channel ipo (z\, z 2 , Z3), the closed 
(atom-molecule) ip c (z,z m ), and the product channels 
Lp d (z,z m ). Elimination of the closed and product chan- 
nels (like in [T2, [I3|) results in a single equation for 



Kid = J2 iTrfVo (°> °> °) + 7««>e (°> °) 
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in terms of the wavefunctions of the open and closed 
channels in the origin 
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Here the transition amplitude X (q,po) is a solution of 
Eqs. (J6]) with Tm (k) HJ), where the energy-dependent 
interaction strength is approximated by 



U eS (E c ) = U a 



%l 



E c + iO - D c 



(12) 



of the two-channel model [l7| ■ 

The coefficients in Eq. (fTTj) are expressed as 
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The same deactivation rate can be obtained within an- 
other model, assuming deactivation due to three-atom 



interactions, 
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FIG. 2: Three-atom correlations for the open |ipo (0, 0, 0) | 2 
(dashed line) and closed \tp c (0, 0) | 2 (solid line) channels, as 
well as the elastic reflection probability P rc f (0 — > 0) (dotted 
line), calculated as functions of the scaled collision energy 
£2 for U a = and D c — 0. The dot-dashed line displays 
lad ¥ a (z) *t ( z ) $t ( z ) q, a ( 2 ) q, d ( z ) |^ (o, 0,0) | 2 within the LLMG model. 



3|£> d | 
4m 



1/4 



*«»(«) *d (*X*%h o ± = [2/i 
pressed as [12} 



)] , the parameters can be ex- 



in place of two-atom ones ©. 

If the non-resonant case (5 = 
is reduced to the LLMG model 
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= 0) the problem 
151 ], which has an 

exact Bethe-ansatz solution. In agreement with [I], Kid 
of Eq. (fTl"]) is proportional to the three-atom correlations 
(fTj), which does not describe suppression of deactivation. 

A physical sense of this effect can be explained by a 
simple 2B analogy. Consider a collision of ID atom and 
molecule with coordinates y and x, respectively. The 
wavefunction of this system has the form of ip (y, x) — 
exp (ipo (y — x)) + i?exp (ipa\y — x\), where R is the re- 
flection amplitude. In the LLMG model atom-diatom 
reflection is forbidden, R — 0, and the 2B analogy leads 
to <p (0, 0) = 1 + R = 1. The 3B model results in Eq. (TJ, 
which has a non-zero value too since the three atoms can 
approach each other. 

However, reflection becomes allowed when integrabil- 
ity of the LLMG model is lifted e.g. by a Feshbach res- 
onance 17J. Moreover, reflection becomes the dominant 
channel at low collision energy, preventing approaching 
of the atoms and leading to suppression of deactivation. 

This hypothesis is confirmed by numerical calculations. 
As the coefficients 7 are independent of the collision en- 
ergy and resonance detuning in both models, ((3]) and 
(fl"4|) . the behavior of the deactivation rate (fTTj) is de- 
termined by the three-atom correlations for the open 
|v?o (0, 0, 0) | 2 and closed |(^ c (0,0)| 2 channels, see Fig. 
[2j The results are expressed in terms of dimension- 
less parameters: the non-resonant interaction strength 
u, the collision energy 62 = 3po/ (4m£>o), and the de- 
tuning b, where Dq is the energy scale. In the quasi-lD 
regime, when p^/m <C lo± and the 3D elastic scatter- 
ing length et3£> does no exceed the transverse waveguide 
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b = D c /D = {(i (B - B 



with (3\ = 1 — Ca3D/a±_ and = 1 + 

C'a 3D fi (B — B — A — w±/m) 2 / (2a ± ujj_A). Here A is 
the phenomenological resonance strength, p, is the differ- 
ence between the magnetic momenta of an atomic pair 
in the open and closed channels, B — B is the detun- 
ing of the external magnetic field B from its resonant 
value B , C & 1.4603 @, and C w 1.3062 [12]. For 
example, in a waveguide with the transverse frequency 
lu± = 50 x 27r KHz, collisions remain quasi-lD for the 
collision energy less then 2ApK. In vicinity of b = 0, Eq. 
(fI5]) gives D = 3pK, u = 0.024, and dB/db = 34 mG for 
the Na resonance at 907 G and Do = 2.5/iK, u = 0.12, 
and dB/db = 27 mG for the 87 Rb resonance at 1007 G. 
Figure [2] demonstrates that the correlations and, there- 
fore, the deactivation rate K\u decrease proportionally 
to the collision energy for slow collisions, when the to- 
tal elastic reflection is approached. It is surprising that 
non-integrability leads to the same low-energy behavior 
of deactivation rate of broad molecules as in collisions of 
structureless particles 13, [Hj]. At rather high collision 
energy, when the elastic reflection probability decreases, 
correlations in the open channel follow to the LLMG 
model. Deactivation suppression persists in a wide range 
of the resonance detunings and the non-resonant interac- 
tion strengths (see Fig. [3]). 

The foregoing results are related to high deactivation 
energies \D&\. The case of low deactivation energies can 
be considered within the two-channel model 17| . A Fes- 



4 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



FIG. 3: Three-atom correlations for the open | <po (0, 0, 0) | 2 
(dashed lines) and closed \ip c (0, 0) | 2 (solid lines) channels, 
calculated as functions of the scaled detuning for £2 = 1 x 10 -4 
and u — 1 (pluses), u = — 1 (crosses), and u = (no symbols). 
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FIG. 4: Deactivation rate coefficient Km (solid lines) and 
elastic reflection probability P rc f (0 — > 0) (dashed lines) , calcu- 
lated as functions of the scaled collision energy €2 for b = 0.4, 
u = 10 (crosses), b = —2, u = —2 (pluses), and b — —4, 
u = — 1 (no symbols), within the two-channel model. 

hbach molecule can have two bound states at U a < 
and D c < 2\g\ 2 /U a , when Eq. (|10p has two real positive 
solutions, Ki > kq. Collision with third atom can lead to 
transitions between the corresponding states, and foreign 
deactivation product states are not more necessary. This 
case is exactly described by Eqs. g]), ©, © and (TP2)) . 
The deactivation rate is given by Eq. ([S]), which includes 
now a singe term (n — 1) only. It is again proportional 
to the collision energy for slow collisions (see Fig. 2]) . De- 
activation suppression correlates with the increase in the 
elastic reflection probability in this model too. 

In summary, the integrable LLMG model does not pre- 
dict suppression of deactivation of broad ID molecules, 
demonstrating their difference from compact molecules. 
The deactivation becomes suppressed when the symme- 
try of LLMG model is lifted by Feshbach resonance. This 
effect is predicted by the two-channel model, as well as by 
two multichannel models with different interactions, ([3|) 



and (fT4| . Dibosonic Feshbach molecules in atom waveg- 
uides become thus relatively stable, like difermionic ones 
in free space. Thus, both the presence [13] and suppres- 
sion of certain processes are among the observable effect 
of non-integrability. 

The author is very grateful to Yehuda Band for stim- 
ulating discussion. 
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